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We calculate the masses of taste non-Goldstone pions and kaons in staggered chiral perturbation 
theory through next-to-leading order in the standard power counting. The results can be used to 
quantitatively understand taste violations in existing lattice data generated with staggered fermions 
and to extract the u, d, and s quark masses and Gasser-Leutwyler parameters from the exper- 
imentally observed spectrum. The expressions for the non-Goldstone masses contain low-energy 
couplings unique to the non-Goldstone sector. With two exceptions these enter as coefficients of 
analytic terms; all the new couplings can be fixed by performing spectrum calculations. We report 
one-loop results for the quenched case and the fully dynamical and partially quenched 1+1+1 and 
2+1 flavor cases in the chiral 5(7(3) and 5(7(2) theories. 
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I. INTRODUCTION 

The masses of the up, down, and strange quarks are 
fundamental parameters of the Standard Model, and the 
low-energy cotiplings (LECs) of chiral perturbation the- 
ory (xPT) [H, parametrize the strong interactions at 
energies soft compared to the scale of chiral symmetry 
breaking, A^. By fitting lattice QCD data to xPT, 
the light quark masses, Gasser-Leutwyler couplings, and 
other LECs can be determined 0, 0] . 

The results of lattice QCD calculations contain dis- 
cretization effects that in principle must be taken into 
account, cither before fitting to xPT or in the fits them- 
selves. For sufficiently small lattice spacings, lattice arti- 
facts perturb the continuum physics d, Q , and the effects 
of these perturbations at energies much less than K-^ can 
be described by an effective field theory, lattice chiral 
perturbation theory 0] . 

Staggered fermions possess an exact chiral symmetry 
at nonzero lattice spacing and are computationally cheap. 
However, in practice discretization effects known as taste 
violations are significant even with Symanzik improve- 
ment. In Rcf. [7[ Lee and Sharpc introduced the xPT 
for a single flavor of staggered fcrmion coupled to SU (3) 
lattice gauge fields. Working to leading order (LO) in a 
dual expansion in the quark masses and lattice spacing, 
they showed that the staggered pion spectrum, includ- 
ing taste violations, respects 5*0(4)7- taste symmetry and 
confirmed this prediction of the xPT by comparing to lat- 
tice data generated by using unimproved and improved 
versions of staggered fermions [1, Q • 

Motivated by unsuccessful attempts to describe lattice 
data by fitting to standard continuum xPT Aubin 
and Bernard generalized the Lee-Sharpe Lagrangian to 
multiple flavors to describe the effects of lattice artifacts, 
including taste violations, in the pseudo-Goldstone bo- 
son (PGB) sector [ll| - [l3j . They used the resulting stag- 



gered chiral perturbation theory (SxPT) to calculate the 
masses and decay constants of taste Goldstone pions and 
kaons (flavor-charged states) through one loop, includ- 
ing the leading chiral logarithms, which enter at ncxt-to- 
leading order (NLO) in the dual expansion [H, [3]. The 
results were used to successfully describe lattice data and 
factored in phenomenologically successful calculations of 
quark masses, meson masses, decay constants, form fac- 
tors, mixing parameters, and other quantities 0, [HI, [lB~ 

In Ref. [1^ Sharpe and Van de Water enumerated the 
complete NLO Lagrangian of SxPT and used it to pre- 
dict relationships between taste breaking splittings in the 
PGB masses, decay constants, and dispersion relations. 
The NLO Lagrangian breaks 50(4)7- to the lattice sym- 
metry group and contributes to the masses of the PGBs 
terms at NLO in the dual expansion. 

Lattice QCD calculations with staggered fermions are 
conducted by taking the fourth root of the fermion deter- 
minant to eliminate remnant doubling from the sea 
The conjectured relationship of the rooted staggered the- 
ory and QCD has implications that can be numerically 
tested. The rooting is systematically incorporated into 
SxPT using the replica method [12, [H, Ho]- Rooting 
leads to unphysical effects at nonzero lattice spacing. Fol- 
lowing the arguments of Refs. (3ll - l43} . we assume that 
the unphysical effects of rooting vanish in the continuum 
limit and that SxPT with the replica method correctly 
describes the effects of rooting that enter soft pion pro- 
cesses at nonzero lattice spacing. 

Here we calculate the masses of the taste non- 
Goldstone pions and kaons through NLO (one loop) in 
SxPT. The results can be used to determine the up, 
down, and strange quark masses, the Gasser-Leutwyler 
couplings, and other quantities by confronting lattice 
data generated with rooted staggered fermions. Con- 
sistency between lattice data and the SxPT description 



2 



of the taste non-Goldstone sector would constitute addi- 
tional numerical evidence for the conjectured relationship 
between the rooted staggered theory and QCD. 

In Sec. |TT]we review the formalism of S^PT, in Sec. IIIII 
we calculate the self-energies of the taste non-Goldstone 
states, Sec. IIVI contains the resulting one-loop corrections 
to the masses, and in Sec. |V] we discuss the results and 
note directions for future work. Appendices El |B1[C1 and 
|D]respectively contain a derivation of the power counting 
through NLO, a discussion of the taste symmetry break- 
ing induced by NLO analytic terms, details of the loop 
calculations, and details of the calculation of the coeffi- 
cients of the chiral logarithms. 



II. STAGGERED CHIRAL PERTURBATION 
THEORY 

Here we briefly review the formulation of S^PT [3, [l^l , 
recalling relevant differences between the staggered the- 
ory and its continuum counterpart [l|. 
we consider the 4-1-4-1-4 theory of Refs. 
metrics and degrees of freedom of S^PT are recalled in 
Sec. Ill A[ the extension of the power counting in Sec. Ill B( 
and the Lagrangian in Sec. Ill CI 

Pedagogical treatments of lattice xPT are given in 
Refs. [4J, 143. Investigations of the foundations of SxPT 
were reported in Refs. js^ Is^ - Isg! ] . 



For simplicity 
il. Thesym- 



A. Group theory and degrees of freedom 

In the continuum limit, the chiral symmetry of the 
A+4:+A theory is SU{12)l^ SU{12)r. Assuming sponta- 
neous breaking to the vector subgroup, the 143 pseudo- 
Goldstonc bosons (PGBs) in the adjoint irrep of SU{12)v 
can be classified according to the continuum flavor-taste 
subgroup 5t/(3)F X SU{4)t: 

SU{l2)v D SU{3)f X SU{4)t (!) 
143^ (8, 15)©(8, 1)©(1, 15) (2) 

Discretization effects break the direct product of the con- 
tinuum chiral symmetry and Euclidean rotations to a di- 
rect product of the lattice chiral symmetry and hypercu- 
bic rotations [l2j: 

U{l)v X SU{12)l X SU{12)r x 50(4) (3) 



^ U{3)l X U{S)r X X 5VF4,diag) 



(4) 



U{3)i X U{3)r is the lattice chiral symmetry of three fla- 
vors of staggered fermions. Its appearance ensures that 
a nonet of the PGBs becomes massless in the chiral limit 
even at nonzero lattice spacing; by definition, these are 
the taste Goldstone states. The U{l)a is not to be con- 
fused with the anomalous axial symmetry of the theory; 
the flavor-taste singlet meson receives a large contribu- 
tion to its mass from the anomaly and is not among the 
PGBs. 



The Clifford group r4 is a subgroup of taste SU{4:)t; 
r4 is generated by the Hermitian, 4x4 matrices 
{Ca^iCi^} = 2<5Ai!y O- SWAAiag is the group of hypercubic 
rotations embedded in the diagonal of the direct product 
of Euclidean rotations and the remnant taste 50(4)t 
that emerges at energies soft compared to the scale of 
chiral symmetry breaking : 



50(4) X SU{A)i 



C 



5W4,diag (5) 

50(4) X 50(4)t (6) 



where Eq. ([6]) represents the symmetry of the staggered 
chiral Lagrangian at leading order. 

The emergence of 50(4)t implies degeneracies among 
the PGBs. The fundamental rep of SU{A)t is a spinor 
under 50(4)t, and the SU{A)t adjoint and singlet of 
Eq. © fall into five irreps of 50 (4) t El: 



5C/(4)t D 50(4)t 

15^ P© A©T( 
1 ^ / 



V 



(7) 
(8) 
(9) 



The 50(4)t irreps are the pseudoscalar, axial vector, 
tensor, vector, and singlet (or scalar), respectively. The 
flavor-nonet taste-pseudoscalar PGBs are the taste Gold- 
stone states. Among them are the pious and kaons of 
Refs. [H, The taste singlet 77' receives a large con- 
tribution to its mass from the anomaly and can be inte- 
grated out of the theory. 

At nonzero quark masses, the continuum chiral sym- 
metry is softly broken to SU{l2)v, SU{9,)v x 5C/(4)y, 
or SU{A)v X 5?7(4)y x 5C/(4)y for three degenerate, two 
degenerate, or three non-degenerate flavors, respectively. 
Noting the anomaly contribution in the taste singlet sec- 
tor and assuming the taste singlet PGBs are degener- 
ate with their physical counterparts in the continuum 
limit, we can use the isospin, strangeness, and contin- 
uum vector symmetries to deduce the degeneracies be- 
tween the remaining (taste non-singlet PGBs) and the 
physical states. 

In doing so we consider the target continuum theory 
with 1-|- 1-1-1 flavors, in which there are 12 valence quarks 
(and 12 ghost quarks), but the fourth root has reduced 
the number of sea quarks from 12 flavors to three. We 
also restrict our attention to PGBs constructed exclu- 
sively of valence quarks. The resulting deductions from 
symmetry represent one of the simplest testable impli- 
cations of the correctness of the rooting conjecture. We 
can also use them to check our S^PT calculation of the 
masses of the taste non-Goldstone pions and kaons. 



B. Power counting 

The standard power counting of SxPT is a straight- 
forward generalization of that in the continuum the- 
ory 0, |46| . The Lagrangian is expanded in a series of 
local interactions perturbing the low-energy theory about 
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the chiral and continuum limits, and observables are cal- 
culated in a dual expansion in the quark masses and lat- 
tice spacing. 

The order of an operator in the Lagrangian corre- 
sponds to the number of derivatives, quark mass fac- 
tors, and lattice spacing factors in the operator. The 
symmetries of staggered fermions ensure the leading lat- 
tice artifacts are i^'(a^), and derivatives always appear 
in pairs Q- Let ri,p2, rim, and be the number of 
derivative pairs, quark mass factors, and (squared) lat- 
tice spacing factors in an interaction. The general form 
of the Lagrangian is 



J^" — ^ ^ -S^2n — -5^2 + -554 -|- . . . 

n=l 

CJO 



(10) 

(11) 



where 



n„2 -I- rim + ria^ and J/f^ 



This organization of the Lagrangian is consistent with 
the expectation that contributions at leading non-trivial 
order will be 



^(pVA^) « mm,/A^) « ffia'Al). 



(12) 



This power counting is appropriate to data generated on 
the MILC asqtad coarse lattices (a « 0.12 fm); on finer 
lattices or with more improved versions of the staggered 
action, the discretization effects are often smaller. 

Feynman graphs are functions of external momenta pi , 
the quark masses mq, and the lattice spacing a^: 



(13) 



where the amplitude ^ is related to the §-matrix as 
follows: 



(14) 



Rescaling pi , ruq , and to smaller values in accord with 
the power counting in Eq. (I12p . we have 



i,mq, a^) — ^{y/epi,emq, ea^), 



which leads to |4j 



^{yfepi^emq^etP') = e^^{pi,mq,a'^), 

oo 



(15) 

(16) 
(17) 



A derivation of Eq. p7p is given in Appendix [Xj is 
the number of loops in the graph, and N2n is the num- 
ber of vertices from operators in ^2n- From Eqs. 
and (fT7)) . we see that loops and diagrams with vertices 
from higher order interactions are suppressed at small 
momenta, quark masses, and lattice spacings. 

For any given observable we first consider all graphs 
with D = 1 (LO), then those with D ^ 2 (NLO), and 



so on. At leading order {D = 1) the only solutions to 
Eq. ((ni) have iV2n = for n > 2 and Nl = 0; i.e., only 
tree graphs with vertices from the LO Lagrangian are 
allowed. At NLO (D = 2), the solutions have N2n = 
for n > 3 and cither A^4 = 1 or Nl = 1; one- loop graphs 
with vertices from the LO Lagrangian and tree graphs 
with at most one vertex from the NLO Lagrangian are 
allowed. In Sec. IIIII we use these observations to write 
down the graphs contributing to the PGB self-energies 
through NLO in the dual expansion. 



C. Lagrangian 

The Lagrangian is constructed of the PGB fields 0, 
quark mass matrix Af , derivatives, and taste matrices 
in accord with the symmetries of the terms in the effective 
continuum Symanzik action 0, fl^ . [29j . The exponential 
parametrization is a convenient way to include the PGBs. 



SU{12)lx SU (12) r: Y.^ LY.R^ (18) 




where L, i? G SU{12)l,r and 

= ^ ® T'' (19) 

K (20) 

SaJ 

e {^5, «^m5, < J^), Cm. 6}. (21) 

The index a runs over the 16 PGB tastes in the 15 and 1 
of SU{A)t, the 0° are Hermitian 3x3 matrices, and the 
are Hermitian 4x4 generators of U{4)t, chosen (up to 
phases) as members of the Clifford algebra generated by 
the matrices With this choice for the T", the S'0(4)t 
quantum numbers of the PGBs are explicit. 

We follow Refs. [H, [2^ in including the flavor-taste 
SU{12)v singlet in the Lagrangian, so S G U{12). An 
additional mass term in the Lagrangian accounts for the 
anomaly contribution to the mass of the SU (12) v singlet. 
Taking this mass correction to infinity at the end of the 
calculation properly decouples the SU{12)y singlet and 
yields the desired results [3 13 ■ 

At leading order in the expansion of the Lagrangian, 
there are three classes of interactions: operators with 
(np2,nm,na2) = (1,0,0), (0, 1,0), and (0,0, 1). We have 

-SfLo = ^Tr(a^I]5^I]t) _ L/2Tr(ME + MI]t) 
o 4 

+ ^{Ui + Di + Sif + a\^ +'^') (22) 



where 



M 



' TO„ 
TOd 

m, 



(23) 
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^7 is the identity matrix in taste space, and the trace (in 
flavor-taste space) is ordinary; we use the rephca method 
of Damgaard and Sphttorfi^ to generahze the results of the 
4+4+4 theory to the partially quenched case [30j . 

The term proportional to TOq is the contribution from 
the anomaly, and the potentials and "i/' break SU (4)^ 
to SO{\)t. They are 

CiTr(ef^S^^"^E^) 
+ C6 5^Tr(^(:)l]e("^)Et) 

+ C3i^[Tr(e^I]e(")l]) + /i.c.] 

V 

+ C4i5][Tr(C(^)sei'^)l]) + /i.c.] (24) 



- = C2yi ^[Tr(e(")l])Tr(^(")S) + h.c\ 

V 

+ C24^[Tr(c(5's)Tr(e(",'l]) + /..c.] 
+ C5yi^[Tr(C(")E)Tr(e)st)] 

+ C5AiE[T^-(^^5^^)Tr(d"'st)] (25) 

V 

where T''^") = T'^'^) ^ i^f^T" in the 4+4+4 theory 
and r° is given in Eq. [2TJ /a is the identity matrix in 
flavor space. The potentials are derived by mapping the 
operators of the mass dimension six effective continuum 
Symanzik action into the operators of xPT. The rem- 
nant taste symmetry of Lee and Sharpe emerges because 
contributions to the potential from 50(4)T-breaking op- 
erators in the Symanzik action are suppressed in the 
low-energy effective flcld theory by powers of the four- 
momenta of the PGBs. The derivation of the potentials 
and the restoration of taste 5'0(4)t symmetry are de- 
scribed in detaU in Refs. 0,[ll[3. 

At NLO, the Lagrangian operators fall into six classes: 
(np2,n„,n,2) = (2,0,0), (0,2,0), (1,1,0), (1,0,1), 
(0, 1, 1), and (0, 0, 2). The flrst three contain terms anal- 
ogous to those in the Gasser-Leutwyler Lagrangian 
The last three contain the terms enumerated by Sharpe 
and Van de Water [2^. The Gasser-Leutwyler terms of 
SxPT that contribute to the PGB masses at NLO are 

^GL- L4Tr(a^I]ta,,E)Tr(xtS + xS^) 
+ L5Tr(a^S]ta^S(x^S + Et^)) 

-L6[Tr(xtl] + xSt)]2 

- L8Tr(x^Sxtl]) + h.c. (26) 

where x = '2fJ.M . 

Many operators in the Sharpe- Van de Water La- 
grangian contribute at NLO, but only a handful break 
the remnant taste 50(4)7' to the hypercubic subgroup 
SW4 of the lattice theory [2^ . We use the symmetries 



of the Sharpe- Van de Water terms to deduce the form 
of their contributions to the masses; as discussed in Ap- 
pendix |Bl the explicit results of Sharpe and Van de Wa- 
ter for the S'0(4)T-breaking contributions to the flavor- 
charged PGB dispersion relations restrict the number of 
independent parameters in these contributions to only 
three. 



III. SELF-ENERGIES OF FLAVOR-CHARGED 
PSEUDO-GOLDSTONE BOSONS 

The symmetries protect the flavor-charged PGBs from 
mixing. For reasons discussed in Sec. IIII Al below, here 
we describe the calculation in the rest frame. In terms 
of the self-energy i;(p|) at p = of the state (j>'^y {x ^ 
y, x,y e {u,d,s}), 

Ml ^ml + E(-A/2), (27) 

where is the tree-level (LO) mass, and is the 
(exact) mass to all orders. Noting that the perturbative 
expansion of £(^4) begins at NLO and expanding £(^4) 
in a Taylor series around p\ ~ gives 

Ml = ml + E(~m2) - E(-A/|)E'(-m2 ) + . . . 
= + E(-m|) + NNLO, 

and the NLO correction to the mass is the leading con- 
tribution to S(— mp. 

In Sec. IIII Al we consider the Fcynman graphs entering 
the expansion of the self-energies at NLO. In Sec. IIII Bj we 
outline the calculation of these graphs, and in Sec. IIII CI 
we present a condensed version of the results for the 
4+4+4 theory. 

A. Diagrammatic expansion 

The power counting of Sec. IIIBI the Lagrangian of 
Sec. Ill CI and the definition of S constrain the diagrams 
entering the NLO mass corrections to three types. Ex- 
panding the LO Lagrangian through and the NLO 
Lagrangian through we write 

(28) 

where cr"^"" corresponds to the sum of connected tadpole 
diagrams (Fig. [T]), a'^'^^'^ corresponds to the sum of dis- 
connected tadpoles (Fig. [5]), and cr"""' corresponds to 
the sum of tree- level diagrams (Fig. |3]). The tree- level 
diagrams are analytic in the quark masses and (squared) 
lattice spacing, while the loops contribute the leading 
chiral logarithms. 

The 4-point vertices in the tadpole graphs are from 
the ^(</)*) terms in the LO Lagrangian of Eq. ([2^ . and 
the 2-point vertices in the tree-level diagrams are from 
the (^{(fP') terms in the NLO Lagrangian of Eq. ([25)1 and 








FIG. 1. At next-to-leading order, tadpole graphs contribute 
the leading chiral logarithms. The vertices are from the lead- 
ing order Lagrangian of Eq. ((22}, and the propagator repre- 
sents the connected part (first term) of Eq. (|29p . 




FIG. 2. Disconnected tadpoles enter in the flavor-neutral, 
taste-singlet, -vector, and -axial channels. The open circle 
represents the second term of Eq. (|29p . 



FIG. 3. At next-to-leading order, tree-level graphs contribute 
terms analytic in the quark masses and (squared) lattice spac- 
ing. The vertices are from the Gasser-Leutwyler and Sharpe- 
Van de Water Lagrangians of Eq. (|26[l and Ref. [2^ . 



readily checked against the taste Goldstone case [12| . 
Sec. nil CI contains a condensed version of these results, 
which can be straightforwardly generalized to the par- 
tially quenched 1-f H-l theory and other cases of interest 
(SecHVl). 

We calculate the loops without extracting the vertices 
by summing over the flavor and taste indices in the ^((/)^) 
terms in the Lagrangian. Instead we combine the expres- 
sions for the (tree-level) propagators of the flavor-charged 
and flavor-neutral PGBs and perform the Wick contrac- 
tions before summing over the vertices. In Ap- 
pendix [C] we detail the calculation of the contributions 
from the mass and a?''W vertices. 



the NLO Lagrangian of Sharpe and Van de Water [29l |. 
The disconnected propagators (in the graphs of Fig. [5]) 
are from quark-level disconnected contributions to the 
tree-level, flavor-neutral propagators in the taste singlet, 
axial, and vector channels |12l |. 

The one-loop graphs break taste SU{A)t to the rem- 
nant taste SO{^)t of Ref. 0, the tree- level graphs 
from the Gasser-Leutwyler Lagrangian respect SU{A)t, 
and the tree-level graphs from the Sharpe- Van de Wa- 
ter Lagrangian break SU{A)t m two stages: Terms of 
^(a^m,) and ^(a") break SU{4)t to S0{4)t, while 
terms of ^{a^p^) break the spacetime-taste symmetry 
SO{4) X 50(4)t to 5W4,diag [II- 

The one-loop graphs respect spacetime SO{4) rota- 
tions, and the corresponding contributions to the self- 
energies, CT^°" and cr'***^, are functions of p^. The 
S'0(4)T-breaking analytic terms of ff{a^p^), however, 
cannot in general be written as functions of p^: The dis- 
persion relations are distorted at nonzero lattice spacing 
by the taste violations. To extract the masses one may 
consider the self-energies in the rest frame. In this case 
the self-energy may be written as a function of pi, the 
square of the temporal component of the momentum. In 
Appendix [B] we recall the form of the 5'0(4)T-breaking 
corrections to the dispersion relations [29j . 



1. Propagators and vertex classes 



Expanding the LO Lagrangian of Eq. through 
^) yields the propagators [H, Hlj. They are 



cab 



SilSjk 



(29) 

where i,j,k,l e {u,d,s\ are flavor SU{3)f indices, a,b 
are taste indices in the adjoint irrcp, and 



Sa 



where 



{q'^ + Ia){q^ + La) 
{q^ + Ua){q^+Da){q' + Sa) 

{q^ + ^?Mq^ + Va){q^ + V'a) 



Si = Aml/3, (5^. = 0, 65=0 



A' 



(30) 



(31) 
(32) 



and the names of mesons denote the squares of their tree- 
level masses. For X G {/, J, i, S}, 



(33) 



B. Calculation in 4+4+4 theory 

For the 4+4+4 theory of Ref. [Hj , we outline the cal- 
culation of the graphs in Figs. [U [2l and [3l After writ- 
ing down the propagators and vertices (Sec. IIII B we 
present results for the loops corresponding to each class 
of vertices (Sec. IIII B 2|) . These intermediate results are 



where nix is the mass of the quark of flavor x e 
{i,j,l,u,d,s}, while for X € {tt", 77, ry'}, the squares of 
the tree-level meson masses are the eigenvalues of the 
matrix 



Ua + 5a 5a 5a 

5a Da + 5a 5a 
5a 5a Sa + 5a , 



(34) 
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The squared tree-level mass of a flavor-charged meson 
(l^ly (x ^ y) is 



Pt = -{Xt + Yt) = Ai(m^ + my) + a^At, 



(35) 



where X e {[/,£», 5} and x ^ y e {u,d,s}. The 

hairpin couplings S'yj^ and taste splittings are com- 
binations of the couplings of the LO Lagrangian [l^l : 



/ 16 

~ -p{C2V — C^v) 

^'a = -T^{C'1A - C^a) 
J 



and 



A/ 

Ay 



16 

7 

16 

7 



2 (4C3 + 4C4) 



2 (Ci + 3C6 + C3 + 3C4) 



16 

7 



2 (4C6 + 2C3 + 2C4) 



Aa = ^(Ci-f 3C6 + 3C3 + C4) 
Ap = 0. 



(36) 
(37) 

(38) 
(39) 

(40) 

(41) 
(42) 



We note the symmetry of Eqs. under simultane- 

ous interchange of vector and axial taste labels {V ^ A) 



and the coefficients C3 and C4. The axial taste matri- 
ces i^^5 generate the Clifford algebra; as a consequence, 
the LO Lagrangian is invariant under simultaneous in- 
terchange of the fields 0^ O- —0'*^ and the coefficients in 
the pairs (C3,C4), {C2v,C2a), and (CsyjCs^)- 

The minus sign arises in the unitary transformation 
connecting the bases {$^} and {iC^s}- ^^^^ ^ such that 



(43) 
(44) 



U is unique up to a phase: 



(45) 



and the other taste generators are invariant under U. 

Noting the diagrammatic expansion and the conserva- 
tion of 5*0(4)7-, we see that taste vector and axial fields 
must always appear in pairs in the calculation of the self- 
energies through NLO, and the minus sign in Eq. ([i^ is 
inconsequential. The results in the taste singlet, tensor, 
and Goldstone (pseudoscalar) channels must be invariant 
under interchange of the coefficients in the above pairs, 
while the results in the taste vector and axial channels 
must switch. This symmetry provides a check at each 
stage of the calculation. 



J 



Expanding the LO Lagrangian of Eqs. ([22]), ((24|). and ([25]) and keeping terms of ^{(p'^) gives 11 classes of vertices. 
From the kinetic energy we have two classes, from the mass terms wc have one, and from the potential a^'^, we have 
four: 



^Tr(9^I]5^S]t) = 
-a2CiTr(e^")se^")s]t) = 

lJ.<v 

-a2C3i^[Tr(^(")SCl")S) + /i.c.] = 



^ Ta,cd {d,cf>t, ^% d,<j^li - d^r,^ d,^% (pl,(pi) 



48/2 
12/4 



{jabcd + 3T5ab5cd ~ 4:T5a5bcd) (f^ti^t^^k'Pkl^'li + 



12/ 



T974 'y^X'^abcd + iT^jab 



12/4 



y^(T"afced + 3x^5 ,ab,vb,cd ~(" '^'^ub,a,ub,bcd) ^ij4'j 



(46) 
(47) 
(48) 
(49) 

(50) 
(51) 



where the indices a, &, c, d run over the 16 tastes in the 15 and 1 of SU{A)t, and Tabc--- = ^■[{T'^T T'^ ■ ■ ■) are traces 
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of products of (Hcrmitian) taste matrices. The four operators of the potential ar^' each give one class: 



T(^2yi[Tr(^(")l])Tr(e)l]) + h.c\ - C5yi[Tr(e(")S)Tr(C(")Et)] 



^Jl{C2V - Cay) ^T, 



2a 



^ E(C2vi[Tr(C(")l])Tr(e^S) + /i.e.] + C5yi[Tr(e(")E)Tr(C(")Et)] 
a? 

= -^{C2V + Cay) ^ T.abT.ed + . . . 

' ^(C2^i[Tr(C^5's)Tr(e<:^E) + /i.e.] + C5Ai[Tr(e^5^S)Tr(e(^^ 



a 



Finally, expanding the Gasser-Leutwyler Lagrangian of Eq. (|26p through ^{(fP') gives 

8i4 



8i 



L4Tr(a^Sta^S)Tr(x^S + x^^) = -^9^0-^9^'/' ^5 + ^5 + ^5) + 
i5Tr(a^Eta^S(x^E + Stx)) 
-L6[Tr(x^E + xSt)]2 
-L8Tr(x^Ex^S) + /i.e. 



f 2 



0" /5 4(C/5 + + 55 



8£8 

p 



(52) 

(53) 

(54) 

(55) 

(56) 
(57) 
(58) 
(59) 



The normalization of Eqs. ([5^ and ([55)) differs from that in the continuum xPT by additional factors of 4. These 
factors are systematically canceled by factors of 1 /4 when using the replica method [l^, [3 1131 to arrive at the results 
for the 1+1+1 flavor theory. 



2. Results by vertex class 



We consider external fields and 0*^, where t is 
the taste index, x ^ y , and we use the renormalization 
scheme of Refs. 0, [lj]- For the tadpole graphs with ki- 
netic energy vertices (Eq. (06])), we find 



r 



ij,kl 



ra 

^ij.kl 



(61) 
(62) 



12/ 



_l_ K°- ~) — 2fl°* K°- 

IX ' yi,iyj xx,yy 



Substituting for the propagators and performing the in- 
tegrals for the connected contributions gives 



1 



12/- 



E 



(60) 



1 



(47r)- 



,(63) 



where i = m, d, s runs over the flavors in the loops, a is where Q runs over the six flavor combinations xi, yi for 
the taste of mesons in the loops, O""^ = jTatab = ±1 if « S {u,d,s}, Qa is the squared tree-level meson mass 
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with flavor Q and taste a, and 

1{X) = x{\xiX/K^ + 5i{\fXL) 



(64) 



for any squared meson mass X. The finite- volume cor- 
rection 5i{y/XL) is [ni 

^i(v^^) = ^L ' (65) 



and 5i{\/ XL) in infinite volume. Here the temporal 
extent of the lattice is assumed infinite, L is the spatial 
extent of the lattice, \fXL is assumed large for all values 
of the quark masses, K\ is a Besscl function of imaginary 
argument, and the momentum in units of 27r/L, n G 1? ^ 
is summed over all integer components except n = 0. 
The diagrams with the mass vertices (Eq. (|47p ) are 



1 



12/- 



■E 



i^xy ^yi)^yi,iy 



(66) 



where — ^{mi + nij). Substituting for the prop- 
agators and performing the integrals for the connected 
diagrams gives 

T^E (i^E(^^ + Q5)^(Qa) 



(67) 



where X5 = m^^ and — ni^y. 

For the graphs with vertices in Eqs. (gS]), (gH]), ([SO)), 
and (|5ip . from the potential we obtain 



1275 4- 



xx,yy 



(68) 



where 



J2Cb{5 + 36'''' e'"^ - 46*''^'' 6*^* - 46'''' e'"^) (69) 



= — ^ Cb{l + W''e''* - 2e^''6''' - 2r^6'''5), 

b^I 



(70) 



and the sum over b includes all tastes appearing in the 
taste matrices of the vertices from ; i.e., all tastes ex- 



cept the taste singlet. The coefficients Ct are the cou- 
plings of the vertices from : 



Ci if 6 = 5 

Ce if 6 e 

C3 i^be{^l} 

Ci if 6 e {^5}. 



(71) 



Substituting for the propagators and performing the in- 
tegrals for the connected contributions gives 



972 



12/ 



(2^)4 



(I^E^w^) 



(72) 



For the contributions from the (hairpin) vertices of 
Eqs. (1221) and we have 



6/ 



^ ^ Sa{2 + e''')Y,{Klxx+Klyy), (73) 



where a runs over the taste vector and taste axial irrcps, 
V = {fj,} and A = {/Lt5}. Proceeding as above, this result 
becomes 



6/ 



1^ 5:^.(2+^*) 



1 



(47r)- 



;il{Xa)+l{Ya)) 



(74) 



where Xa and Ya are given by Eq. ((33|). The connected 
and disconnected pieces of this result can be combined 
into a single disconnected piece by using the identity [12 1 



1 



q^ + Ia 



Sa 



(75) 



where la and La arc given by Eq. (|33p . ? G {u,c?, s}, 
and a e {y, A, /}. This result follows immediately from 
the form of Dfj obtained by treating the ^(0^) terms of 
the LO Lagrangian that are proportional to Sa as vertices 
and summing the resulting geometric series for the flavor- 
neutral propagators: 



Sa 



1 



(q^+IaW + Ja) 



1 



(76) 



q^ + La 



The equivalence of this form and that given in Eq. ((30)) 
was demonstrated in Ref. [l^ for general partially 
quenched theories. With the use of the identity Eq. (|75p . 
Eq. ([71)1 becomes 

- ^ E^(2 + n / ('^' + (77) 
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For the loops from the remaining two vertices of 
Eqs. dSll) and we find 



with flavor Q and taste a, and Pa is given by Eq. (|35l) 
Setting = —Pt = — P5 — a^At, we have 



32/ 



-Ef E -.Wa.*(l + r^)V™' (78) a-(-P,) = -^E(«'^-"«'(^* + ^<'))^(Q" 

a \beV,A I a,Q^ 



where 



^ 16 I + Csy if6eM ,7q^ 



Performing the integrals gives 



(80) 



For the tree- level graphs with the vertices in Eqs. (|56p. 
ST]). (fSSj) . and ([Ml), we have 



i|(2L6P5 + i4p')4(C/5 + i55 + 55) (81) 



16 

-^i2LsP5 + L5P^)P5 



(82) 



while the tree-level graphs with vertices from the Sharpe- 
Van de Water Lagrangian may be parametrized by intro- 
ducing LECs corresponding to the irreps of S'0(4)t and 
SW4: 



Ifi 

—a^£/tP5 + 4(C/5 + D5 + S5) 
f 



(83) 



where the coefficients £/t, iMt, and are degenerate 
within the SO{4:)t irreps, and the coefficients 'ift are 
degenerate within the SWi irreps. The symmetry of 
the calculation under interchange of the (valence) flavors 
X <r^ y implies the &(a?mq) terms are proportional to P5 
or the sum C/s + + 6*5 . 



C. Results in 4+4+4 theory 



The results in Eqs. f63 |) . (|67 )) . ([721). ([77 | . and ([80 1) are 
the one-loop contributions to the expansion of the (nega- 
tive of the) self-energies of the flavor-charged PGBs of 
taste t S {/,^,/xzy(/x < :^),/i5,5} in the 4+4+4 the- 
ory. Collecting the connected contributions and factoring 
-l/(47r/)2 gives 



a.Q 



32 ^ 



E ^6^afatra&t(l+0"') 



l{Pa). (84) 



where Q runs over the six flavor combinations xi, yi for 
i € {u,d,s\, Qa is the squared tree-level meson mass 



32 ^ 



E ^bTabtTabtil + 0'''') 



.b<£V,A 



l{Pa)- (85) 



For the Goldstonc case, i = 5 and Eq. (|69|) with Eqs. ([SS 
through (|42| imply 



while 



Aa5 = Aa, 



T-a65rab5(l+0 =Oif 5e 



(86) 
(87) 



so the connected contributions vanish identically [l2| . 

The chiral logarithms are degenerate within the 
50(4)7- irreps; summing over the values of a within each 
irrep, we rewrite Eq. (j85p : 



a = —a 



A con 

Y.VbfI{Pb)^^Y.^{Qb)Y (88) 
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where 



^ ^ E E ^bTaUTau{\ + Q'''') (89) 



aeB 



Ag-^^ ^(A„t-(At + A„) 



(90) 



The indices B and F designate the S'0(4)t irreps, B,F E 
{I,V,T, A, P}, t £ F, and we conveniently abuse the 
notation to define 



Xb = Xa for a e B and X e {P, Q}, 



(91) 



which is possible because the taste splittings are degen- 
erate within irreps of SO{A)t- 

The coefficients S'j^p and A'^p are linear combinations 
of the couplings in the potentials and , respectively. 
Equivalently, the coefficients A'^p are linear combina- 
tions of the taste splittings A^ , and the coefficients Sgp 
are linear combinations of the couplings ujv a (defined in 
Eq. dZD). 

Explicit results for 5'^p and A'^p are given in Tables |T] 



PF 



^PF 



0. A 



PF 



follows 



and ini We note that S 
from the symmetry of the summand of Eq. (|90|1 under 
interchange of the indices a, t and the vanishing of the 
connected contributions in the Goldstone case, Ag'p — 0. 
The symmetry a t does not exist in the summand of 
Eq. ([89l) . so the relation (5pJ? — appears non-trivial. In 
general the symmetry a O i of the summand of Eq. ((90|) 
implies 



(92) 
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TABLE I. The coefficients 5%°^ defined in Eq. dHlJ. To obtain 
(5s°p, multiply the entry in row B and column F by 16//^. 



B\F 


V 


A 


V 








A 








T 


3(C2A + Csa) 


3(C2i'' + C^v) 


P 








I 


C2V + C^v 


C2A + CsA 



B\F 


T 


P 


/ 


V 


2{C2A + Cba) 





4(C2v + Csv) 


A 


2{C2v + Csv) 





4(C2A + Csa) 


T 











P 











I 












where Np^^^-f is the dimension of irrep F{B). Eq. is 
useful for checking the results in Table HT] Eq. ([M]) , with 
the coefficients in Tables U and is our final result for 
the connected tadpoles in the 4+4+4 theory. 

Collecting the disconnected pieces from the one-loop 
results in Eqs. tmd ^ gives 



12 



■E 



(27r)4 



a2A,i-a2A,j(Z?^, + D^V 
-26'"y + (2(1 - e"*') + p'^*) g2 



+ (2(i + r*) + p'^*)p5 

+ 2a2A;, + (2 + p'^*)a2A, 



where 



I -4(2 + 0-*) ifa^/ 
^ [0 ifa = /. 

Setting = — Pf = — P5 — a^Af, we find 

(47r)2^ f d^q 



(93) 



(94) 



12 



■E 



(27r)4 



a'Aa 



a^(A,+A,) (7^S. + Z?^%) 



(2(1 - ei''*) + p"*) g2 ^ 1^2(1 + 26*"*) + p'^*) F5 
2a2A:., + (2^°* A* + (2 + p"*) A, 



(95) 



For the Goldstone case, i = 5 and Eqs. (|69| . ((TO]) , and 
imply 



A' 



a5 



(96) 



and Eq. (|95|) reduces to the result of Ref . [11] . 

The sum over a receives nonzero contributions from 
the 50(4)7' vector, axial, and singlet irreps. Summing 
over a within each irrep, we can write 



disc 



(47r)2 r d^q 



12 



(2^)^ 



2( -I2P5 - Qq''>^VF + a'A*|f 



+ {V^A)+a'Ar^HDi, + Dly) 
+ 2(3P5 + a^Af^^)Diy 



(97) 



where 



Af |f =Y,(Kt + O'^'At + (1 + p'^V2)A„ (98) 

aeB ^ 



aeB 



for t E F. i^BF is the number of taste matrices for irrep 
B E {V,A} that commute with the taste matrix cor- 
responding to t E F. The values of i>bf a-re given in 
Table HVl 



The coefficients A 



disc 
BF ' 



like the coefficients A' 



BF^ 



are 



linear combinations of the taste splittings Aq. In Ap- 
pendix |D] we show that 



Ag|f = forP 



A disc 
VI 

A disc 
AI 



A^7 - 24A 
AX/" - 24A 



L,V,A 

V 
A- 



(100) 
(101) 
(102) 



The latter two relations follow from the identity 

a;, - A,i-2(A, + Ai) ifr*-l. (103) 



Eqs. (jlOip and (|102p provide non-trivial checks of the 
results for Aff in Table |TTT1 Eq. with the coefla- 

cients in Tables HIl IHIi and llVi is our final result for the 
disconnected tadpoles in the 4+4+4 theory. 

Taking into account the hairpin couplings, taste split- 
tings, and coefficients (5|°|5, Ag'^, and Af|f in Tables [fl 
im and mil we see that the loops, the results in Eqs. 
and (|97p. are invariant under the symmetry 



C3 O C4 
C2V ^ C2A 



(104) 
(105) 
(106) 



or, more brieffy, under 

V ^ A (107) 
in accord with the observations following Eq. (HSI) above. 
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TABLE II. The coefficients Ag'|5 defined in Eq. To obtain A^i^^, multiply the entry in row B and column F by 96/ f. 



B\F 


V 


A 


T 


P 




V 


4Ci +C3 + 9Ci + 6C6 


4Ci + 3C3 + 3C4 + 6C6 


2C3 + 6C4 + 86*6 





4C3 + 12C4 


A 


4Ci + 3C3 + 3C4 + 6C6 


4Ci + 9C3 + C4 + 6C6 


6C3 + 2C4 + 8C6 





12C3 + 4C4 


T 


3C3 + 9C4 + 12C6 


9C3 + 3C4 + 12C6 


6C3 + 6C4 + 16C6 





12C3 + 12C4 


P 

















I 


C3 + 3C4 


3C3 + Ci 


2C3 + 2C4 





4C3 + 4C4 



TABLE III. The coefficients Aj^"/ defined in Eq. ((98|. To 
obtain A*!?, multiply the entry in row B and column F by 
96/ f. 



B\F 


V 


A 


V 


-3Ci - 6C4 - 


3C6 


-Ci - 6C4 - 


9Ce 


A 


-Ci - 6C3 - 


9Cg 


-3Ci - 6C3 - 


3C6 


I 


C3 + 3C4 




3C3 + Ci 







T 


P 


I 


V 


-2Ci - 8C4 - eCe 





-4Ci - 12C6 


A 


-2Ci - 8C3 - 6Ce 





-4Ci - 12C6 


I 


2Ci + 2C4 





4C3 + 4C4 



TABLE IV. The numbers vbf of taste matrices for irrep B G 
{V, A} that commute with any given taste matrix for irrep F. 
Vbf appears in row B and column F. 



B\F 


V 


A 


T 


P 


/ 


V 


1 


3 


2 





4 


A 


3 


1 


2 





4 



Collecting the analytic contributions to the self- 
energies from Eqs. ([81]), dMl), and ([831) gives 

a'^""'(p2) ^ i|(2Lg P5 + Lipl) 4(C/5 + -D5 + 55) 

+ ^(22.8^5+ is P4) ^5 

+ ^a2(^tP5 + 4(C/5 + + 

+ '^tpl + ^ta^). (108) 
Setting pI = ~Pt — —P5 — a'^At, we have 

^anal ^ ^ ^^^^ _ ^^)p^ ^(f;^ + i?5 + ^5) 
Ifi 

+ -{2Ls~L,)Pi 
16 

+ -j^a^c^tPa + 4(C/5 + + ^5) 

+ ^ta^), (109) 



where we have absorbed terms proportional to a^At into 
the coefficients (ft, J^t, and ^#t. The first two lines of 
Eq. (|109p correspond to the continuum result and are the 
same for all tastes. In the last two lines, the coefficients 
are degenerate within irreps of 50(4)^, while the co- 
efficients (ot and are degenerate within irreps of SW4. 
The exact chiral symmetry implies that ^5 = = 0. 
Setting t = 5 in Eq. (jlOOp then yields the result of 
Ref . [I4] . In Appendix [B] we recall the results for the 
dispersion relations of Sharpe and Van de Water [2^ ; the 
S'0(4)T-breaking contributions to the S't and terms in 
Eq. (|109p come from only three operators in the Sharpe- 
Van de Water Lagrangian. The consequences for fitting 
strategy are discussed in Sec. [V] 

Eqs. (IMl), and ((1(7^1) are useful starting points for 
deriving results in various cases of interest. In Sec. II VI we 
use them to deduce results for fully dynamical, partially 
quenched, and quenched theories. 



IV. RESULTS 

The results in Eqs. and pTO| must be mod- 

ified to account for (partial) quenching [5^, [U and the 
fourth root of the staggered fermion determinant [3| • The 
replica method of Ref. [s^ allows us to generalize to the 
partially quenched case. The replica method also allows 
us to systematically take into account the fourth root of 
the staggered determinant [13,, i32„ ,34 ' .36.] . 

The effects of partial quenching and rooting in 
Eqs. (IMl), (EZl), and pH^]) are easily summarized: The 
valence quark masses and ruy are no longer degener- 
ate with the sea quark masses m„, rud, and rris, a factor 
of 1/4 is introduced in the second term of Eq. ((88|) . the 
eigenvalues of the mass matrix in Eq. (|34p are replaced 
with the eigenvalues of 

Ua+6a/4. 6a/4. 5a/4 \ 

6a/4: Da+5a/4: 5a/i , (110) 

5a/4: da/4: Sa+Sa/y 

and terms in Eq. (jl09p that are proportional to the sum 
of the sea quark masses are multiplied by 1/4. 

Accounting for quen ching the sea quarks is equally 
straightforward [ll|. [T2|. [soj : The second term of Eq. §^ 
is dropped, and the disconnected part of the propagator. 
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in Eq. (l30l) , is everywhere replaced with 

^quench 



D 



a, quench 



where 



ucnch ^ I 4(^0 + aq^)/^ if a = / 

if a ^ /. 



(Ill) 



(112) 



The one-loop contributions to the self-energies in the 
fully dynamical l-l-l-l-l and 2-1-1 flavor cases in the chiral 
SU{3) theory are in Sec. IIV ATI One-loop contributions 
for partially quenched 1+1+1 and 2+1 flavor cases of in- 
terest are in Sec. IIVA21 In Sees. IIVBII and IIVB2I we 
write down the analogous (fully dynamical and partially 
quenched) one-loop contributions in the chiral SU (2) the- 
ory. Sec. IIV A 31 contains one-loop contributions for the 
quenched case. 



A. SU{3) chiral perturbation theory 

1. Fully dynamical case 

Introducing a factor of 1/4 in the second term of 
Eq. ((88)1 . we obtain the connected loop contributions in 
the fully dynamical 1+1+1 flavor case, 



a = —a 



\C07l 

E('^bf^(^b) + ^E'(Q^)1' (113) 



48 



where Pb is the LO squared mass of a flavor-charged 
PGB (pl^y with xy valence (anti)quarks {x y), and Qb 
is the LO squared mass of a PGB (pf^ with z£ valence 
(anti)quarks, where z G {x, y}, £ G {u, d, s}, and the sum 
over Q runs over the six flavor combinations formed by 
pairing the possibilities for z with those for £. Setting 
xy = ud, us, ds gives the results for the 7r+, K^, and 
We have 



_con 2 \ ^ / ccon 

cr^+ =-a ObfI(t^b) 



A con 
^BF 



1{Ub) 



+ 2/(4) + 1{K+) + 1{Db) + l{Kl)) ^ (114) 

(A con , 
5bfKKI) + ^(i{Ub) 

l{n+) + 21{K+) + liK%) + liSB))^ (115) 

. / Aeon , 



„con J2 

CT^o = —a 



+ 



1{Db) + iml) + 1{K%) + 1{Sb))\^ , (116) 



where the squared tree-level masses of the flavor-charged 
mesons are 



^{■niu + nid) + a^As 
fi{mu + m 



Kb = ^J-{md + rus) + a^A^. 



(117) 
(118) 
(119) 



The integrals of the disconnected pieces, in Eg. (|97L 
can be performed by using the residues of Ref. [I4I to 
expand the integrands. For the 7r+ we find 



disc 



1 

12 



^VF^'v E 

X 



+ R^Dl°nn'{Xv)l{Xv) 



2'Y^i-12Ti+ + QXvvvF + a^^ 



X 



disc 
VF 



a^5'yR^.JXv)l{Xv) + {V^A) 



Va-E 

X 



R§^o^iXj)liXj) + i?g^o,(X,)/(X,)^ 



'-Arj^^)j2R'.o,{Xj)iiXj) 



X 



(120) 



where we decoupled the flavor-taste singlet, rj'j, in the 
taste singlet channel by taking tiiq 00 before expand- 
ing the integrands. 

In Eq. (jl20p we introduce a few naively perverse but 
extremely convenient abuses of notation. First, in each 
sum over X, the residue in the summand determines the 
values of the index X . For example, in 



T 

X 



(121) 



the index X takes the values in the set {!7, 7r°, 77, 77'}. 
When the summation over X is factored from sums of 
residues specifying different sets of values for the index X , 
as in the first and second lines of Eq. (|120p , we first dis- 
tribute the summation symbol and then use the residues 
to specify the values of X in each sum. Second, the 
SO{4)t irrep specified in the argument of a given residue 
applies to all masses appearing in the residue. For exam- 
ple. 



R 



DS 

UiT^riri' 



(Uv) 



{Dy - Uv){Sv ^ Uy) 
{n°^-Uy){vv~Uy){v'y~Uy)-' 



(122) 



where we continue denoting squared tree-level masses by 
the names of the mesons. In general the residues are 



R 



BiB2-B„ 



(123) 



where X e {^1,^2, 
50(4)t irrep. 



IlBi=^xi^iF - Xf) ' 
,B„} and F e {V,A,I} is the 
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The results for the K+ and K° may be obtained by 
permuting U,D,S in the residues of Eq. (|120|) and re- 
placing TT"*" with K~^,K'^. They are 
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Sv{Ru.o^r,'iXv)l{Xv) 



{Xv)l{Xv) 



a2,5i,i?^o„„,(Xy)/(Xy) + (F ^ A) 



(124) 



-(3X+ + a^Ar/)i?-o,(X,)Z(X,) 



and 



„disc 



-E 
12 ^ 



- 2 (^-12X0 + QXyvvF + a'^vF^ x 
a25;.i?^„ + {V^A) 



2 A con / tdUS 



-a'A 



IF 



K'.oJXi)l{Xi) + R''sX{Xi)l{Xi) 



+ -(3^0 + a^ATF^R^iXiyiXj) 



(125) 



In Eqs. ((nni), (dH]), and ([T^ . the masses tt^, 77s, and 
are the eigenvalues of the mass matrix in Eq. (|110p 
with rriu rud- 

For physical values of the quark masses, strong isospin 
breaking is a small correction to electromagnetic isospin 
breaking, and 2+1 flavor simulations have proven very 
useful . Setting xy — ud, us and rriu — rrid in Eq. (jll3p 
gives the connected contributions for tt, K: 



performing the integrals after setting to„ = m^. We find 



12 



X 



2^( -127r5 + 6Xy 



- a'{A^y-^ + Aff))aH'yR%^,{Xv)l{Xv) 



+ iV^A) 

+ ^{3n, + 2a'Arpn(lK^i)~lKvi) 



, (128) 



disc 



12 



X 



a*A^y"-6'yJ2(R%^,iXv)liXv) 



- Rs^^,{Xv)l{Xv 

- 2 (-12X5 + &XvvvF + a^AyF 

X ^ 

aH'yRr,r,'{Xv)l{Xv) + {V ^ A) 



IF \ 2^M + KSi) 



UiK,~\a'ArF]Km) 



(129) 



where ttb, rjB, and rj'^ {B e {y, A, I}) are the eigenvalues 
of the mass matrix in Eq. (jllOp with m^j = rrid, and we 
used the relations of the tree-level masses in the taste 
singlet channel to simplify the associated residues: 

R^^^im)^^ Rl_^(^j) = -^ (130) 
RUSi)^3 RUvi) = -'2- (131) 



In the continuum limit, only the taste singlet contribu- 
tions to the disconnected loops survive. Taking the con- 
tinuum limits of Eqs. ([T^ through ([T^ . we recover the 
one- loop results of Gasser and Leutwyler 



A con 
^BF 

24 



Aeon , 

a" ( 5%°^ 1{Kb) + [21{^b) + 51{Kb) 

(127) 



2/(7rs) + 1{Kb) 



(126) 



48 
1{Sb) 



In the 2-1-1 flavor case, the disconnected contributions 
are most easily obtained by returning to Eq. (j97l) and 



quenched case 



The connected contributions in the partially quenched 
1-1- 1-1-1 flavor case have the same form as the connected 
contributions in the fully dynamical 1+1+1 flavor case, 
Eq. pi3p . The difference is that the valence and sea 
quark masses are, in general, non-degenerate: rnxjiny ^ 
{mu,md,ms}. 

For the disconnected contributions in the 1+1+1 flavor 
case, keeping all quark masses in Eq. (j97l) distinct and 
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performing the integrals as before, we find 
1 



_disc 



a AypSyyRx^Orir,'{Xv)liXv) 

+ i?^fo^,„, iYv)l{Yv) + J2 {dI%.^,^x {Zv)l{Zv 
z 

+ D^yX^,^y{Zy)l{Zy))^ 

+ 2 ^ (^-12P5 + &ZyVyF + A^'ff ^ X 



a'5'yR'i'il,,^^,{Zy)l{Zy) + {V ^ A) 

4 
— ( 
3 



Y,(Dl''J^^x{Zj)l{Zj)+D^%^^y{Zj)l{Zj 



(132) 



where 
D 



and 



(Xf) (133) 



1{X) = -[\nX/K^ + 1] +5:i{^/XL). (134) 



The finite- volume correction 5^{\/XL) is [11 1 
5^{^/XL) = 2Y,KQ{\n\^L), 

n=iO 



(135) 



and S3{\/XL) — > in infinite volume. 

A non-trivial special case of Eq. (|132l) occurs for 



my. We have 



disc 

x=y - ^2 



2 1^-12X5 + a2(A^°^ + A^'|f ) j X 
a 5y ( R^^o^^, 

{Xy)l{Xy) + {Zv)l{Zv. 

12a2j(,l.y;^(^i?5[^^;^,(Xy)(;(Xv) - Xyl{Xy)) 

(Zy);(Zy) 
Z ^ 

+ [V^A) + ^(3X5 + 2a2A™") feo^ rx,)[(X,) 



^i?^^cf,,;,(Z,)?(Z,) 



(136) 



The masses ir^.rjB, and 77^ (S e {V', A, /}) in Eqs. (fT32|) 
and (jl36p are the eigenvalues of the mass matrix in 

Eq. (inni). 



To obtain the connected contributions in the partially 
quenched 2-1-1 flavor case, we set to„ = in Eq. (|113p . 
To obtain the disconnected contributions, we set viu — 
nid in Eq. (j97p and consider the two cases rrix 7^ rriy and 
rrix — my separately. We find 



disc 



1 

12 



a'A'y"^6'y(Rf^^,{Xy)l{Xy) 



Rf^^,{Yy)l{Yy) + Y.iD\%^,^x{Zy)l{Zy) 



{Zv)l{Zv. 



2 ^ (^-12P5 + QZyvyp + a'A^'l?^ : 



^'Wy,,'(^i^)^(^v) + (V^^A) 



+ -a'Af^ \R%{Xi)l{Xi) + i??-^,(17)Z(17) 
+ Y.{^x%,xiZinZi) + i?^^,y(Z,);(Z,) 



(137) 



and 



^disc 



1 

12 



2 (^-12X5 + a"(A^°;? + A^*|f ) I X 

a^5'y { R\%^,{Xy)l{Xy) + D%^,^x{Zy)l{Zy) 
-12 a^^V^VF {Rl%^,{Xy)(l{Xy) - Xyl{Xy) 

~ X! ZvD''xr^^, x{Zv)l{Zv 
z 

+ {V^A) + \{iX^ + 2a2AJ-) ( Rf^{Xi)l{Xi) 



Y.D%,x{Zi)l{Zi\ 



(138) 



The masses tt^, 773, and rj'g [B e {V,A,/}) appearing 
in Eqs. (|137p and (|138p are the eigenvalues of the mass 
matrix in Eq. (|110l) with m^ — ma- 



3. Quenched case 
The connected loop contributions are 

a'°^ ^-o'Y.^bfKPb). (139) 

B 

To obtain the disconnected contributions, we consider 
Eq. ([W)) with the replacement — > where 
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^a.quench ^j^^^ ^T^-^T^. We have 



—disc 



12 



+ 12 a Oyi^vF — 



Xv-Yv 

+ {V^A) 



(^{ml - aXi)l{Xi) + aZ(X/) 



+ + 3(3P5 + a'A™")x 

(m§ - ayj)^(yj) - {ml - aXi)l{Xi) 
X,-Y, 



(140) 



and 



cr. 



disc 



1 

12 



12 a 



+ (F ^ ^) + -(3X5 + 2a"Ar/)x 



({ml - aXi)l{Xi) + al{Xi)^ 



(141) 



where in Eq. (|140p we substituted for the residues, 



Rxy{Xb) — —Rxy{Yb) 



YR-Xf 



(142) 



The loop contributions to the pion and kaon masses in 
the case of three non-degenerate quarks and in the isospin 
hmit can be obtained from Eqs. (jl39l) . (|140p . and (|14ip 
by appropriately choosing m^; and my. 



B. 5(7(2) chiral perturbation theory 

Expansions about the SU{2) chiral limit are often 
better behaved than expansions about the SU{3) chiral 
limit. The correspondin g y PT was developed in Ref. [l| 
and extended in Refs. |52l . [53| to describe results ob- 
tained with rooted staggered fermions. To date SU{2) 
SxPT analyses have been restricted to the taste Gold- 
stone sector [ssj - lssj . Beginning with the loop contri- 
butions from SU (3) SxPT, we write down correspond- 
ing SU{2) SxPT loop contributions to the taste non- 
Goldstone PGB masses. 



1. Fully dynamical case 
From Eqs. pT4)) . (fTT5| . and (|TT6| . we have 

(A con , 
SBFl{^t) + ^(l{UB) + 2l{^%) 

+ 1{Db))^ (143) 

. . \con , , 

(145) 



'K+ 



B 



^con 



A con y 



B 



To obtain the disconnected loop contributions in the 
1-1-1+1 flavor case, we consider Eqs. ([TMl) . 
and Eq. (HZl). We find 



disc 

^17 + 



1 

12 



X 



a^ATF5'vY.[R^^o,{Xv)l{X^ 



RDTTOrj{Xv)l{Xv) 



2 \ disc 



2 I -127r+ + GXyiyvF + A 

X ^ 

a^6'yR^o^{Xv)l{Xv) + {V^A) 
2a2A— (r§.o{Xi)1{Xj) + RlAXi)l{Xi) 



X 



+ 4(3715+ + a^Af/)Z«) 



(146) 



^disc 



12 



a^A'y°^S'yJ2RS.o,{Xv)l{Xv) 

X 

-12^ a'S'yR%^{Xy)l{Xy) + {V^A) 

X 

+ 2a'ArF-Y^uAXi)l{Xi) 



X 



(147) 



^dis( 



1 

12 



a'A'y-^S'yY Rd.o,{Xv)1{Xv) 

X 

-12J2 a'5'yR%{Xy)l{Xy) + {V^A) 

X 

+ 2a'Ar^Y.^DAXi)l{Xi) 



X 



+ 6(!7z -<)/(<) 



(148) 
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The taste vector and axial residues in Eqs. (I146p - (jl48p 
can be simplified using the tree-level masses: 



-isgn5's)y{DB-UBy + \ia^S'^r (149) 

+ {sgnS's)y{DB-UBy + \ia^S'sr (150) 
ior B = V,A. We have 



Ru-nOniUB) = 


4 


RuTrO^iA) = 


2 


RuTrOi^iVB) = 


2 


RDiT<>r,iDB) = 


4 




2 


RD7,'>rjiVB) = 


2 



(1 + sin/3s) 
(1 - sin /3s) 



-(1 + sin /3s) 

Id 



- COS /3s 



2 

a 

2 

R^o^iA) = ^(1 + sin/^s - COS 
1 
2 



R'^o-qiVB) = i (1 - sin/3B + cos 



/Sb) 
/3s), 



where 

sin /3b = (sgn(5^) 



^{DB~UB? + \{aH'B? 

cos /3b = (sgn(5g)^ ^ ^ 

v/(i^B-C/B)2 + i(a2^i,)2 

In the isospin limit, /3b = 0. 

The connected loops in the 2+1 flavor case are 



(151) 

(152) 

(153) 

(154) 

(155) 

(156) 

(157) 

(158) 

(159) 
(160) 

(161) 
(162) 



'K 



rcon 
'BF 



A con 
'^BF 

12 



Kt^b) 



A con 
^BF 
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^(7^5), 



(163) 
(164) 



and the disconnected loops are 



1 

12 



4(-127r5 + a2(A^^.°^ + A*|f) 

{liny) - liT^v)) 

+ 2Avvf{'^vK'^v) -VvKvv)) 
+ iV^A) 



+ 4(3^5 + 2a2AJ].")/(^/) 



(165) 



^disc 

'"^ =12 



2a^A^°^(/(^y) - l{7^v)) - I2a'5'yl{r^v) 



+ {V ^A) + 2a2A™"/(7r/) 



(166) 



The masses in Eqs. pl5| . (ITM| . and ([TBS)) are 



■KB = 2/iTO„ + a" /S.B V -B 

2 A a'^i 
775 = 2^TO„ + a As H — 



(167) 

Be{V,A}. (168) 



All mesons circulating in loops in the SU (2) chiral theory 
are pious. 



2. Partially quenched case 



We obtain the connected contributions in the 1+1+1 
flavor case by dropping terms in Eq. (jll3p corresponding 
to loops with a strange sea quark; i.e., the sum over Q 
excludes the xs and ys mesons, and we treat the x and 
y quarks as lightQ To obtain the disconnected contribu- 
tions in the 1+1+1 flavor case, we consider Eqs. (|132p . 



^ Another case of interest would be that of a single heavy (strange) 
valence quark, my ~ rus 2> ~ m„ ^. 
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((T361) . and We have 



_disc 



12 



+ R^^o^{Yv)l{Yv) + Y,(dI%^^{Zv)1{Z^ 
+ D'^%^y{Zv)l{Zv)) 

+ 2 ^ ( -I2P5 + GZvi^vF + a^A'^'pA x 
a'S'yR'^'}^o,{Zv)l{Zv) + {V^A) 

z 

+ 4(3P5 + a^Ar/) ^XY.o {Zi)KZi) 



(169) 



and 



disc 



12 



2 (^-12X5 + a2(A^°^ + A*|f ) j X 

^ z 
l2a'5'yVvF {r"x%{Xv){i{Xv) - ^^[(Xv)) 

T.^y^x%AZv)i{Zv) 



+ {V^A)+ 4(3X5 + 2a2Ar/) i?^^o(X,)Z(X,) 



(170) 



Setting = m^^, we have the disconnected contribu- 
tions in the 2+1 flavor case: 



disc 
x^y,u— 



1 

12 



a'A'y°^6'y (R-x,iXv)liXv) 

R%{Yv)1{Yv)+J2{d3c,,x{Zv)1{Zv) 
z 

D^,^^yiZv)l{Zv))^ 
2 J2 ( -12^5 + GZviyvF + a^A^'ff 



a25^i?^^,,(Zy);(Zv) + {V^A) 
2a^Ar^((l{Xj) + {nj - Xj)l{Xj) 



+ (liYi) + im - I7)f(y/)) 
+ 4(3P5 + a^Aj'P) J2 R\y{Zi)1{Zi) 



(171) 



and 



disc 
^ x—y,u—d 



1 

12 



2 1^-12X5 + a2(A^°^ + A*|f ) j X 
a^S'y(Rl^{Xv)liXv) + J2Dl^,xiZv)KZv) 



rx^{Xv){l{Xv) - XvUXv)) 



12 a^SyVvF [R 
Y.ZvD\^^x{Zv)l{Zv)] + {V ^ A) 

4(3X5 + 2a'AT,^)(l{Xi) + {^i - Xi)l{Xi)^ 
V. CONCLUSION 



(172) 



Our final results for the masses of the flavor-charged 
PGBs through NLO in SxPT are given by adding 
Eq. (|28)) evaluated on-shell to the tree-level (LO) result 
of Eq. ([55]) . These results and others of interest can be 
obtained from those in the 4-1-4+4 flavor theory given in 
Eqs. (IHHl), (EZl), and of Sec. (ImC]) . Applying the 

replica method to reduce the number of tastes per flavor 
from four to one gives the connected tadpole, discon- 
nected tadpole, and NLO (analytic) tree-level contribu- 
tions to the on-shell self-energies. In Sec. IIVI we write 
down the connected and disconnected tadpoles in the 
1-l-H-l flavor and 2-f 1 flavor cases in SU{i) and SU{2) 
XPT. 

For the fully dynamical case with three non-degenerate 
quarks, the results in the SU{Z) chiral theory are in 
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Eqs. (I114p through (I125p . The corresponding results in 
the isospin hmit are in Eqs. (jl26|) through (|129|) . Expan- 
sions about the SU{2) chiral hmit are given in Eqs. (|143|) 
through (ITI51) and Eqs. dM]) through ([TMl) . For the 
quenched case, the resuhs are in Eqs. (|139p through 
()14ip . where the LECs are the quenched counterparts 
of those in the theories with dynamical quarks. 

For the partially quenched case, the connected contri- 
butions have the same form as those in the fully dynami- 
cal case, Eq. (|113p . For three non-degenerate sea quarks, 
the disconnected contributions in the SU{3) chiral theory 
are in Eqs. (|132l) and (|136l) . Taking the isospin limit in 
the sea, the corresponding results are in Eqs. (|137p and 
()138p . The expansions about the SU{2) chiral limit are 
given in Eqs. pMll - pTOl) and Eqs. (IT7T |) - p7^ . 

The LO contributions to the masses break taste 
SU{A)t to taste SO{A)t 0. At NLO the (tadpole) loops 
respect taste SO{4)t, the tree-level counterterms from 
the Gasser-Leutwyler Lagrangian respect taste SU{A)t, 
and tree-level counterterms from the Sharpe-Van de Wa- 
ter Lagrangian break spacetime-taste SO{A) x SU{4)t to 
the lattice symmetry, r4 x S'W4_diag- 

The pattern of taste symmetry breaking is illustrated 
in Fig. m Regarded as functions of the valence masses, 
the LO masses receive corrections at NLO to their slopes 
and intercepts. The chiral logarithms contribute to both 
types of corrections but do not lift degeneracies within 
taste SO{4:)t irreps. A small subset of the Sharpe-Van de 
Water counterterms breaks the S'0(4)t symmetry. With 
HYP-smeared staggered valence quarks on MILC coarse 
lattices, the corrections to the intercepts are smaller than 
the statistical uncertainties jUjH^]; Fig.|3]represents this 
case. The exact chiral symmetry at nonzero lattice spac- 
ing ensures corrections to the intercept of the taste Gold- 
stone (P) mesons vanish. 




mx + my 

FIG. 4. Pattern of taste symmetry breaking in squared PGB 
masses. The PGBs fall into 8 lattice irreps, and the masses 
receive corrections to the slopes and intercepts that lift the 
degeneracy of the taste 50(4)t irreps. The corrections to the 
intercepts are of i^(a*) and very small on typical lattices. 



As discussed in Appendix iBl the S'0(4)T-breaking con- 
tributions to the masses of the flavor-charged PGBs arise 
from only three operators in the Sharpe-Van de Water 
Lagrangian ^2^. As a direct consequence, the S'0(4)t- 
breaking corrections to the slopes and intercepts of Fig. |4] 
depend on only three LECs. One could obtain them from 
the splittings in the slopes of the taste vector, axial, and 
tensor irreps. 

The LECs of the LO taste-breaking potentials enter the 
results in four specific ways: The couplings of the poten- 
tial , Eq. (|24|) , enter only via the tree- level mass split- 
tings and the coefficients of the chiral logarithms from 
disconnected tadpoles; for each taste channel, these co- 
efficients are given in Tables and IIIII The couplings of 
the potential , Eq. enter only via the hairpin co- 
efficients of the Goldstone sector and the two linear com- 
binations of Eq. (f79| . The former arise in disconnected 
propagators, while the latter multiply connected tadpoles 
with valence-valence mesons in the loop; for each taste 
channel, they are given in Table HI 

From Eqs. ([351) -(gH) and Tables HI] and [ml we observe 
that the coefficients of the chiral logarithms from the 
disconnected tadpoles are completely determined by the 
tree-level mass splittings. The tree-level mass splittings 
also determine the coefficients of chiral logarithms from 
connected tadpoles with sea quarks in the loop. Hav- 
ing determined the iS'0(4)T-breaking terms and the LO 
masses, one could perform fits to partially quenched data 
to extract the remaining (two) coefficients of the con- 
nected contributions and the coefficients of the 5*0(4)^- 
preserving analytic corrections at NLO. We note in pass- 
ing that the coefficients of the quenched and SU{2) chiral 
theories are different from the coefficients of the SU (3) 
chiral theory. 

The calculation here can be extended to mixed action 
xPT and to other quantities of phenomenological inter- 
est such as decay constants, form factors, and mixing 
parameters. For example, one could consider HISQ or 
HYP-smeared staggered on asqtad staggered simulations; 
the SxPT for both cases is the same [5a-[6l|. A calcu- 
lation for Bfc is given in Ref. [27[. We plan to calculate 
in the near future the one-loop corrections to the mass 
spectrum of pions and kaons in the mixed action case. 
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Appendix A: Power counting formula 



lattice spacing, 



In Sec. Ill Bl we recalled the standard power counting 
of SxPT. Here we derive the power counting formula, 
Eq. (IT71) . The derivation is based closely on the discussion 
for the continuum case in Ref. [^^. 

We begin by noting that the effective continuum 
Symanzik action contains no operators of mass dimen- 
sion five or seven 12911: 



Ss 



YM 



5*4 



a^Se 



(Al) 



For the following discussion, we assume without proof 
that no operators of odd mass dimension appear at higher 
orders in the Symanzik action. At the end of the deriva- 
tion, we consider the restrictions this assumption places 
on the vahdity of Eq. (fT7| . 

Mapping the operators of the Symanzik action into the 
Lagrangian of S^PT and using it to compute an arbi- 
trary amplitude, we note that dependence on the lattice 
spacing enters via the vertices and the (tree-level) prop- 
agators. By the assumption of the previous paragraph, 
all vertices and propagators depend analytically on a^. 
The Symanzik and SxPT actions are translation invari- 
ant, so momentum conservation holds in SxPT, as in the 
continuum theory. 

The lattice-spacing dependence of the propagators can 
be deduced from Eqs. dSO]), ([32]), ([33]), §5^, (pi)) , 

and Eqs. ([55)1 through (|42l) . The propagators receive cor- 
rections proportional to the hairpin couplings of Eqs. p6p 
and ([57]) and the taste splittings of Eqs. ([551) through 
(|42|) . The former enter only the disconnected parts of 
the flavor-neutral propagators; the latter are corrections 
to the tree-level masses of all the PGBs. 

In any given amplitude, internal lines contribute fac- 
tors of 



(27r)4 
(2^ 



ill) 



(A2) 



SiiSjk- 



+ h{Ia + Ja) 



(A3) 



and vertices contribute momentum-conserving delta 
functions and couplings of the form 



P 



m, 
a 



ea 



the factors from internal lines scale as 



(27r) 



(27r) 



(A5) 
(A6) 
(A7) 



(A8) 



where we changed the variable of integration from q to 
^/eq. Under the same change of variable and the rescaling 
of Eq. (lASp . the momentum-conserving delta functions 
scale as 



8\p + q)^5\^e{p + q)) 
= \6\p + q), 



(A9) 
(AlO) 



while the vertex contribution of Eq. (|A4[) scales as 



e V, 



(All) 
(A12) 



where n = ripi -\- + • Recalling the definitions of 
np2, rim, and and the organization of the Lagrangian 
in Eq. PTI)) . we conclude that Eq. (IA12p implies that a 
vertex contribution scales with a factor of e" if and only 
if the vertex is from an interaction in the Lagrangian 
of Eq. (fTO|. 

A given Feynman graph ^{pi, uriq, a^) has Ni internal 
lines and Ny vertices, where 



E 



No. 



(A13) 



i.e., the number of vertices in the graph is the sum of 
the number of vertices N2n from each term in the 
Lagrangian. Because the total momentum flowing into 
the diagram equals the total momentum flowing out, one 
momentum-conserving delta function does not contribute 
an independent constraint and is factored out of the am- 
plitude. Multiplying the rescaling factors of the inter- 
nal lines, vertex contributions, and remaining Ny — 1 
momentum-conserving delta functions gives 



Ji({pi,mq,a^) 



-'^{pi,mq,a^) 



(A14) 



(A4) where 



where 2np2 is the number of derivatives in the interaction, 
p is an external momentum, q is an internal momentum, 
Hq is the number of internal lines contracted with the 
vertex, rim is the number of quark-mass factors from the 
vertex, is a light (u, d, or s) quark mass, 271^2 is the 
number of lattice-spacing factors from the vertex, and a 
is the lattice spacing. 

Rescaling the external momenta, quark masses, and 



D^Ni- 2{Nv - 1) + 



nN; 



2n- 



(A15) 



The number of loops in a diagram is the number of in- 
dependent integrations after imposing the Ny — 1 con- 
straints from the momentum-conserving delta functions: 



Nl = Ni- {Nv - 1). 



(A16) 
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Using this relation to eliminate Nj from Eq. (|A15|) gives 
the desired result: 



D = Nl-Nv + 1 + Y^ nN2n (A17) 

ri=l 

oo 

= l + AfL + ^(n-l)iV2„. (A18) 



We now reflect on the validity of this result. We as- 
sumed that no operators of odd mass dimension appear 
at any order in the Symanzik action. This assumption is 
known to be true only through mass dimension eight. In 
principle it could be violated at mass dimension nine: 



'SYM 



= 54 + 0^5*6 + 58 + 0^59 + . . 



(A19) 



In this case the rescaling is the same for the internal 
lines and the momentum-conserving delta functions, but 
the possible vertex contributions are different: 



(A20) 



and the expansion of the SxPT Lagrangian can be writ- 
ten 



n=2A,5,... 



(A21) 



We allow for the number of derivatives Up in an operator 
to be odd because they are the only objects in the chiral 
Lagrangian with indices that can contract with those of 
taste matrices to construct operators with an odd number 
of taste spurious. 

The vertex factors now rescale as 



where n = Up + 2n„j + Ua, and we have 



The modified power counting formula is 



D = l + Nr 



E 



n=2,4,5, 



1 Uv„ 



(A22) 



(A23) 



(A24) 



and writing out the solutions to this equation for D = 
1, 2, |, 3 yields the same solutions to the power counting 
relation for D ~ 1,2,3 as before. For D = 1,2, the 
new operators in ^5 do not contribute, and the power 
counting of Eq. (fTT]) is justified through NLO. 



Appendix B: NLO analytic corrections and taste 
symmetry breaking 

The NLO analytic corrections to the PGB masses 
break spacetime-taste 50(4) x 50 (4) ^ to the diagonal 
hypercubic subgroup 5W4.diag of the lattice theory. Here 



we consider the responsible operators and note the pat- 
tern of symmetry breaking in the mass spectrum. 

Sharpe and Van de Water enumerated the NLO La- 
grangian giving rise to the NLO analytic corrections [29j. 
Although many operators contribute at NLO, the vast 
majority respect the remnant taste symmetry 50(4)t; 
only three operators are responsible for the symmetry 
breaking in the masses of the flavor-charged PGBs, and 
all are of type {np2 ,nm,na2) — (1,0,1). For example, the 
operator 



J2 Tr(9^0^^9^0C^ 



(Bl) 



yields the correction to the self-energy for (pl.^ 

J2p,P^^0^\ (B2) 



P 



36V 



which in the dispersion relations breaks (spatial) rota- 
tion invariance and lifts the 50(4)t degeneracies of the 
masses of the PGBs. 

The symmetry breaking corrections to the dispersion 
relations were calculated by Sharpe and Van de Water. 
They have the form 





= f + 


M](l + Kl) 




(B3) 


El 


= f{l 


+ K4 — Ki) + 


My(l + «4) 


(B4) 


Ef 


= plil 


+ — K4) + 




■n) (B5) 


Efj 


= [pI ^ 


-Pj)(l + l^ij 


- + pI + M^{1 


+ Kij ) 

(B6) 


Eji 


^pi + 




- Kii - Kij) + M|.(1 


+ Ku) 

(B7) 


El 


-pIU 


+ Ki5 - K45) 




K,5) (B8) 


El 




+ K45 — Kis) 


+ Afl(l + K45) 


(B9) 


El 


= f + 


Ml{l + K5), 




(BIO) 



where Mt is the flavor-charged PGB mass through NLO 
including all but the taste 50(4)7- violating contribu- 
tions, and we use Kt to denote the St of Ref. [23] ■ Writing 
out the various contributions to Mt more explicitly, 



a' At 



at, 



(Bll) 



where the flrst two terms are the LO result, it is the sum 
of all loop corrections for taste t, and at is the sum of all 
NLO analytic corrections respecting taste 50(4)t. The 
results of Sharpe and Van de Water then imply that the 
NLO 50(4)T-breaking corrections are 



{mly + a'^At)Kt = mlyKt + AtHt. 



(B12) 



Since the corrections are proportional to a? , the NLO 
masses of the flavor-charged PGBs, considered as func- 
tions of TTLx + TUy , receive corrections to their slopes and 
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intercepts. However, these corrections are completely de- 
termined by only three a priori unknown LECs in the 
Sharpe-Van de Water Lagrangian. 

Recalling the 50(4)T-breaking corrections in 
Eq. ([1091), 



P 



(B13) 



we see that Sj contains the 16 a priori independent coef- 
ficients rot and , and that the part of 5t which breaks 
SO{A)t is fixed by only three LECs in the Sharpe-Van 
de Water Lagrangian. The remaining part, which is 
S'0(4)T-symmetric, can be determined by fitting as dis- 
cussed in Sec. El 

In summary, the S'0(4)T-breaking corrections to the 
dispersion relations imply the presence of ff{o?'mq) and 
ff{a'^) analytic corrections to the masses of the corre- 
sponding (flavor-charged) PGBs. These corrections are 
determined by only three LECs, which may be taken to 
be the splittings of the SO{4)t vector, tensor, and axial 
irreps. 



Appendix C: Example loop calculations 



In Sec, nil we wrote down the results for each class 
of vertices in the graphs contributing to the PGB self- 
energies at NLO. Here we detail the calculation of these 
results in two cases: for the (mass) vertices of Eq. (|47|. 
which yield Eq. and the {a^'^) vertices of Eqs. PSI) 
through (inil), which yield Eq. with Eqs. and 

dZOl). 

Perhaps the approach taken here can be extended with- 
out too much difficulty to calculations beyond NLO. The 
integrals associated with two-loop contributions will dif- 
fer from those entering tadpole graphs, but in principle 
one can construct them from the propagator (Eq. (|29p ) 
and the generic forms of loops with (for example) ^(0®) 
vertex classes from the LO Lagrangian. 



Mass vertices 



We begin by considering the vertices in Eq. (jTTl) : 



(CI) 



The corresponding tadpole graphs are the sum of all com- 
plete contractions with the external fields and (/)^^. 
Contractions between a field in the vertex and an external 
field vanish unless the flavor and taste indices match, in 
which case they contribute factors of unity to amputated 
diagrams: 



(C2) 



The remaining fields contract together to give the prop- 
agator of Eq. in the loop, and integrating over the 
loop momentum gives the integral of Eq. (|6ip : 



1 



q^ + hila + Ja] 



There are six ways to contract the vertex fields and ex- 
ternal fields while maintaining the order of the latter, or 
equivalently, while maintaining the order of cc and y. Sup- 
pressing the common factors of the coupling /i/(48/^), 
taste factor Tabcd, and mass m^, we have 



rta,th,cd j^c I cta,tc,l 
xj,iy,yk^jx kl.li ' xj,iy^ 



+ 5. 



ta,td,bc 



.K 



+ 5\ 



■tb.tc.ad 



,kx jk,li 



xj,iy,yi,lx ^^jk,kl ^xk,jy,yl,kx ij,H 

I ctb.td.ac i ctc,td,ab T^€i 

' xk,jy,yi,lx ij,kl ~^ "xl,ky,yi,lx ij,jk 

where 6 is simply a product of Kronecker deltas: 

gab,cd.ef 5,,Skl6mnSr,nS''''6'"^6''f. 

Restoring the taste factor Tabcd and quark mass rui and 
summing over the flavor and taste indices gives 

TttCC my K^.^^^y + Ttbtb my Kl^^yy 

+ Ttbbt my K^i ,^ + Tatta mi Ky^ i 



(C3) 



(C4) 



■ Tatat my K^^ yy ~\- Taatt my Ky^ ^y 



(C5) 



where we used the symmetry of the propagator under 
interchange of the flelds and relabeled dummy flavor in- 
dices. The taste matrices all commute or anticom- 
mutc with one another. Defining 6*°^ such that 

j,aj,b ^ ^ab rj^brj^a y , 



and noting 



{T-f=^i V a. 



(C6) 
(C7) 



Eq. (jCSp becomes 



{2m^ + niy + mi)K^i^^^ 
+ {2my + m.j;+ mi)Kyi,_y 
+29^\m,+my)K^.^^^ 



(C8) 



Restoring the coupling gives the desired result, Eq. (|66l) . 

Contributions from loops with kinetic energy, mass, 
and hairpin vertices (respectively Eqs. (|551). (1571) . and 
(|771) ) cancel against the term with 6'°* in Eq. (|S^ . and 
evaluations of the sums over 6'°* (within irreps of SO{A)t) 
are not needed to arrive at the coefficient of the mass 
term in the result, Eq. (j^Tl) . In general such cancella- 
tions do not occur. The needed sums are collected in 
Appendix |D1 
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Potential a vertices 



Appendix D: Coefficient calculations 



We begin by noting that the flavor structure of the 
vertices from a?"^ is identical to the flavor structure of 
the mass vertices (Eq. (|C1|) ): 



"12/4 

6 



{jahcd + 3r5Qf,5cd - 4T5a56c<i) (t>1j'i>)k^kl4>'ii (C9) 



(CIO) 



12/4 



{Tabcd + iT^abucd + ir^aubcd) <t)1i4>']k<f>ll4'tt (CH) 



+ 4T^5,a,jy5,6c(i) 4>'ij<t> jk^kl4>li- 

(C12) 

Only the overall normalizations (couplings) and taste 
factors differ from those in Eq. (|Cip . Therefore the 
same contractions and corresponding products of Kro- 
necker deltas and integrals that appear in Eq. (jCSp enter 
the calculation, and we can obtain the loops for each 
of the vertex classes in Eqs. (|C9p through (IC12I) from 
Eqs. (|C3|) and (jCSp by taking — )• 1 and replacing Tabcd 
in Eq. (ICSp with the appropriate linear combinations of 
traces from Eqs. dCO]) through (|C12p . 

Noting that the taste factors in Eqs. (jC9p through 
(|C12p all have the same form, viz. 



where 



Tabcd ~t~ '^Tgabscd 4^ Tgagjjcd 



5 for vertices cx Ci 
/it/ for vertices cx Cg 
for vertices cx C3 
z^S for vertices cx C4, 



(C13) 



(C14) 



and recalling Eqs. (IC6p and (jC7p . we find the taste factors 
in Eqs. through (|CT2|) are 



Here we detail one way to calculate the coefficients in 
Tables Hini andHIl The coefficients defined in Eqs. §9^ 
and ([70]) can be computed similarly. As a by-product 
of our calculation, we obtain the intermediate results of 
Tables |V] and I VII which we calculated twice: Once us- 
ing the (anti)commutation relations of the generators to 
count the positive and negative terms in each sum, and 
once using the relation 6°''' = Tabab/4: and a computer. We 
have also checked Table |V] (IVip implicitly, using an inde- 
pendent accounting scheme to arrive at the coefficients 
Ag'^ (AffP) given in Table HI] (|lll| . 

The coefficients S'^gp are defined in Eq. ([89t : 



Sbf-^T. E ^brabtrabt{l + 9^') 



a£B beV.A 



and we note 



|o if {r«,e^(5)} = o. 



so that 



and 



(5^°P = V F 



dip 



UA a F = A 

otherwise. 



(Dl) 
(D2) 
(D3) 



(D4) 

(D5) 

(D6) 
(D7) 

(D8) 



Tttcc + ^T^ttscc - 4.e^'Tststcc = 16(1 - 



46i^V,,_^„ 
46''"(1 + 36''"*r* - 46i^"r*) 
Ttbbt + iTstbsbt - 4.9'''T,tsbbt = 4(1 + 39'''e'' - 405-0^*) 



Ttbtb + STstbstb — "iO-'Tstsbtb 
dbtz-i I onbsnst 



Tatta + STsatsta " i9''T,astta = 4(1 + 3^^^^* - A0''0'-) 
'^atat ~l" '^'^satsat 4^ '^sastat 

= 4r*(l + W'"'9'^ - AO^'e"") 

Taatt + iT.aastt " ^9''' T^asatt = 16(1 " ^^^0^°). 

Replacing the taste factors in Eq. (jCSp with these (and 
taking rrii 1 there) leads directly to Eq. (|55)) . 



Similarly, for B — V we have 



jrcon 
Oyp 



^5]K(r.^0'(l + ^^^) 

+ C^A(T.,^5,t)'(l + 0"'^')] 



4wy if F = / 



<5^°^ = <j 2u;a if F = T 
otherwise. 



(D9) 
(DIO) 
(Dll) 

(D12) 



23 



^AF = :^EK(^^'5.pO'(l + ^''''^) (D13) 



and for B = we have 
1 

32 

+ ^A(r.5,M5,t)'(l + ^?"''^')] (D14) 
lj2[u;y{r,^^,^,f{l + e^'-n (D15) 
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2wy if F = T 



(DIG) 
(D17) 



otherwise. 
Finahy, for B = T we have 

^TF - ^ E [^v(rpA,pt)'(l + 0"'-^ (D18) 



^^,p<\ 



+ C^A(TpA,M5,t)'(l + ^''^'^')] (D19) 
= 17^ E K(TpA.pO'+^^A(TpA,M5.t)'] (D20) 



3wv if = ^ 
= <( 3t^A if F = y 
otherwise. 



(D21) 



These results are straightforwardly obtained by sub- 
stituting for t and counting the number of nonzero 
terms in the sums. For evaluating the traces, the 
(anti)commutation of the generators T", the fact that 
(T")2 = ^j, the orthogonality relations Tr(r°T'') = 4(5°^ 
and the traces over products of Euclidean gamma matri- 
ces are useful. 



TABLE V. Sums for evaluating the coefficients Ag"^. All 
the sums required for the coefficients A^]J can be obtained 
by repeated use of the results in the first three lines. 



t G 


P 


I 


V 


A 


T 




-4 


4 


-2 


2 







-4 


4 


2 


-2 







6 


6 








-2 



















^p,p- 


8 


-8 


4 


-4 





ai^paii5,t 

'^p,p- 


8 


-8 


-4 


4 







-12 


-12 








4 



















Y^ np-jpX ap5,t 


















The coefficients of Table |TT] are defined in Eq. (j90p : 



aeB 



Substituting for the coefficients Aat and taste splittings 
using Eqs. (|69l) and (f86| gives 

24 

A™^ = — ECfc(l + - - fi'^^fi'^^) (D23) 



Writing out the sum over each irrep B gives 



A^S^ = 
48 



(D24) 



+ E '^6(1 + fi'^'''"^'^'''* - 6'^'^'* - 6'^''''') (D25) 
-f-EC3(l + 6i^''6i''* + r*+r") (D26) 

-1- E '^4(1 + 6i'''^'"6"'5^* -I- ei'^'^'* -I- 6*^5^") j . (D27) 



(D28) 



^Tf = ^2 E(^3(l + ^''*) + ^4(1 + 0^''*)) (D29) 
24 



(D30) 



+ E ^6(1 + r'^'^'r'''* - r'^'* - 0""'") (D31) 

-hEC3(1 + ^^''^^' + ^^' + ^'"') (D32) 

-I- E C'4(l - 6i^''6i^^'* 6"^^^* - 0^'') j (D33) 
24 / 

^™F = -«E (D34) 
p ^ 

+ E + 6"^"^'* - r'^^* - 6i^''''') (D35) 
p<i^ 

-I- E C'3(l - + - ^'''') (D36) 

-I- E C'4(l + 6i^''6'^^'* -f 6"^^'* -I- 0^'') j (D37) 
24 / 

^TF^pE E^ex (D38) 

p<A \p<!^ 

+E'^3(i +^^'''^^''* +^^* +^^'''^) (^40) 

E ^4(1 + e^'p^df"^'' + r^'* r-"^) j . (d4i) 
p / 



Inspecting Eqs. (|D29p through (jD41[) . we note the req- 
uisite sums. Their values are given in Table |Vl Using 
Table El to evaluate Eqs. (|D29|) through ((D4T]) yields the 



24 



results of Table |TT1 For F = P {t = 5), all the coefficients 
explicitly vanish, as they must. 



TABLE VI. Sums (in addition to those in Table E]) for eval- 
uating the coefHcients A'^p^. The last three lines can be ob- 
tained from the first two, since ec^'* = 6^*6''\ 



t e 


P 


/ 


V 


A 


T 










12 


-12 





QptQp,pQp.t 

L^p,p. 


-8 


-8 


4 


4 


8 


aP^:t niiv^pX r\pu,t 


-12 


-12 
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QptnupnuS.f 


-8 


-8 


-4 


-4 


8 


Qp5,tgiJ.v,pQpu,t 








-12 


12 





()P^J,tgiJ.pgiJ.S,t 

^P.n- 


-8 


-8 


4 


4 


8 



The coefficients of Table |TTI] are defined in Eq. ([98 



^BF = \ ^'at + ^"'^t + (1 + p"72)Aa ) ■ (D42) 

aeB 



Substituting for the coefficients A^^, taste splittings, and 
coefficients using Eqs. dTO]), ^ and jM]), and jMl) 
gives, for B I, 



24 



nat ( nah nbt 



a€B 

b^I 



and for B = I, 
24 

^ b^I 



(D43) 



(D44) 



For B ^ I, adding and subtracting Eqs. (|D43|) and (|D23p 
gives 



A disc 
^BI 


A con 
— ^BI ~ 


67VbAb (B 7^ /) 


(D45) 


A rfisc 


A con 

— LAprp 




(D46) 


A disc 


A con 
— ^TP ~ 


36At 


(D47) 


A C/ZSC 

^PV 


A con 
— ~^PV 




(D48) 


A disc 


A con 
— -^PA 




(D49) 


A disc 
^VP 


A con 
— ~^VP 




(D50) 


A c^isc 

^AP 


A con 
— ~^AP) 




(D51) 



while for _B = /, comparing Eqs. (jD44[) and (jD23[) gives 

(D52) 



A disc A con 
IF — ^IF ■ 



Eq. (|D45]) implies Eqs. (ITOT|) and they and 

Eqs. (|D50p and (ID51I) can be used to cross-check the re- 
suhs in Table Uni Eq. ((D52]) and Table III give the coeffi- 
cients Afp'^, while the coefficients Ap|f and A^p^ do not 
appear in Eq. (|97l) ; the remaining coefficients in Eq. ((97)) 



are A^*|f and Af|f . 

Writing out the sum over b in Eq. (jD43[) gives 

24 / 

^ aeB ^ 

(D53) 

+ E C6(-l + ^'(6*'^''."^''.* - 61^'''*) -t- r''-'^) (D54) 
+ E C3 (-1 + r*(6i^°6l^* -I- 0^*) - 6"^") (D55) 

(D56) 

and writing out the sums over the vector and axial irreps 
{B = V and B ^ A) gives 
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^VF = « H 2C7i(-l - e^'e^*) 



(D57) 



C6(-l + 6i''*(r'''''6l^'''* - r''-*) -f- e^'^'") (D58) 
Y Cz{-1 + 6"'*(6""''6"^* -I- 61^*) - 6'^'') (D59) 

^ C4(-I - 6»''*(6"^''r5'* - r^'*) -I- 0^'')^ , (D60) 



p 



(D61) 



^ C6(-l + eP^'\9'"'-P9'"'-^ - 6'^'''*) + e^''''') 

(D62) 

^C3(-l-6"'5^*(r''r* + (D63) 

^C4(-l + 0P5,t(^0t^P0t^5,t _^ g,M5,t) _ gMP) 



(D64) 



Examining Eqs. (jD57|) through (|D64p . we note the sums 
beyond those in Table IVl that are needed to evaluate the 
coefficients Ay|f and A*|?. The values of these sums are 
given in Table IVII Using the sums in Tables |V] and I VII in 
Eqs. (|D57| thro ugh (lD64l) yields the result s in T able lllll 
From Eq. p44p and Eqs. (|D53P through (|D56p . we see 
that Ajp = A'^p = Afp" = 0, as necessary for the re- 
sult in Eq. ([97)) to reduce properly in the taste Goldstone 
case. 
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